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Abstract 

Four double-parton scattering processes are examined at the Fermilab Tevatron 
energy. With optimized kinematical cuts and realistic parton level simulation for both 
signals and backgrounds, we find large samples of four-jet and three-jet -|-one-photon 
events with signal to background ratio being 20%-30%, and much cleaner signals from 
two-jet -l-two-photon and two-jet+e'*'e~ final states. The last channel may provide the 
first unambiguous observation of multiple parton interactions, even with the existing 
data sample accumulated by the Tevatron collider experiments. 



There are good reasons to believe that multiple partonic interactions, where two or more 
pairs of partons scatter off each other, occur in many, or even most, pp collisions at the 
Tevatron {y/s = 1.8 TeV). On the theoretical side, multiple partonic interactions are an 
integral part of the eikonalized minijet model |^, which attempts to describe the observed 
increase of the total pp cross section with energy in terms of the rapidly growing cross section 
for the production of (mini)jets with transverse momentum pt > Pr,min — 2 GeV. Sjostrand 
and van Zijl also pointed out that including multiple interactions in the PYTHIA event 
generator greatly improves the description of the "underlying event" in pp collisions. A 
similar result was found recently by the HI collaboration in a study of 7p collisions. 

However, hadronic event generators have many ingredients. This makes it difficult to 
draw unambiguous conclusions from such studies. It is therefore desirable to search for more 
direct evidence for multiple partonic interactions, using final states that are amenable to a 
perturbative treatment. Clearly the cross section will be largest if only strong interactions 
are involved. The simplest signal of this kind is the production of four high— p-p jets in 
independent partonic scatters within the same pp collision |^ (4 — 4 reactions). Since energy 
and momentum are assumed to be conserved independently in each partonic collision, the 
signal for a 4 — »• 4 reaction is two pairs of jets with the members of each pair having equal and 
opposite transverse momentum. Various hadron collider experiments have searched for this 
signature. The AFS collaboration at the CERN ISR reported p a strong signal. However, 
the exact matrix elements for the QCD background 2 — 4 processes were not used and the 
size of the signal claimed was considerably larger than expected. The UA2 collaboration at 
the CERN SppS collider saw a hint of a signal, but preferred to only quote an upper bound 
[^. More recently, the CDF collaboration at the Fermilab Tevatron found evidence at the 
2.5cr level that 4 A processes contribute about 5% to the production of four jets with 
Pt > 25 GeV [0. 

While final states consisting only of jets offer large cross sections, they suffer from severe 
backgrounds. There are three possible ways to group four jets into two pairs. Further, the 
experimental error on the energy of jets with pt — 20 GeV is quite large. Hence even four- 
jet events that result from 2-^4 background processes often contain two pairs of jets with 
transverse momenta that are equal and opposite within the experimental errors. The study 
of "cleaner" final states has therefore been advocated: The production of two pairs of leptons 
(double Drell-Yan production) has been studied in refs. [Q, the production of two J/ if) mesons 
in refs.p, and the production oi a.W boson and a pair of jets in refs.[|l^]. However, in our 
opinion none of these processes is ideally suited for studying multiple partonic interactions. 
Double Drell-Yan production offers a very clean final state, but the cross section at Tevatron 
energies is very small once simple acceptance cuts have been applied. The cross section for 
double J/ifj production is quite uncertain, since it depends on several poorly known hadronic 
matrix elements Finally, VT+jets events can only be identified if the W boson decays 



leptonically, which makes it impossible to fully reconstruct the final state. 
Here we study mixed strong and electroweak final states: 

• three jets and an isolated photon (jjj'y); 

• two jets and two isolated photons (jj'yy); 

• two jets and an e~^e~ pair (jjee). 
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For comparison, we also include 

• four-jet final states (denoted by 4-jet). 

We try to be as close to experiment as possible within a parton level calculation. To this 
end we not only apply acceptance cuts, but also allow for finite energy resolution, and try 
to model transverse momentum "kicks" due to initial and final state radiation. We find 
that the jjj'j final state offers an only slightly better signal to background ratio than the 
4-jet final state does; note that the combinatorial background is the same in these two cases. 
This combinatorical background does not exist for the jj'yy and jjee final states, which offer 
much better signal to noise ratios, at the price of small cross sections. 

The calculation of our signal cross sections is based on the standard assumption |^, ||, 
1^, |l^ that the two partonic interactions occur independently of each other. The cross section 
for a 4 — >■ 4 process is then simply proportional to the square of the 2 — > 2 cross section: 

a(4 ^ 4) = [a(2 ^ 2)] Vao. (1) 

This assumption cannot be entirely correct, since energy-momentum conservation restricts 
the available range of Bjorken— a; values of the second interaction, depending on the x values 
of the first one. We include this (small) effect using the prescription of ref.[^. In the 
eikonalized minijet model [Q ctq is related to the transverse distribution of partons in the 
proton. Unfortunately total cross section data do not allow to determine this quantity very 
precisely. We find values between about 20 and 60 mb, depending on the choice of the 
numerous free parameters of the model. The recent CDF study ^ found do = 24.2t^Q;g mb, 
within the range that can be accommodated in minijet models. We will take ctq = 30 mb in 
our numerical analysis; the results can be scaled trivially to other values of ctq. 

The for us relevant 2 — 2 cross section can be written as a sum of different terms: 



o-(2 ^ 2) = a{pp jjX) + a{pp j-fX) + a{pp 77X) + a{pp e+e X), (2) 



where j stands for a high— p-r jet. Inserting eq.(0) into eq.(|l]) gives a 4 — > 4 cross section 
that sums over many different states; it should be obvious which terms in the sum are of 
relevance to us. Note that this procedure gives an extra factor of 2 in the cross section 
for the production of final states made up from two different 2^2 reactions (e.g, jjj'j) 
compared to those produced from two identical reactions. Partly for this reason we only 
consider jjjj configurations where the two jets are produced in one partonic scatter and 
the two photons in another. The other possible configuration (77^7), where each jet pairs 
up with one photon, also suffers from larger backgrounds, since there are two ways to form 
such pairs. We use leading order matrix elements in eq.(|D, but we include the contribution 
from gg 77, which enhances the total pp 77X cross section by about 50% at y/s = 1.8 
TeV. We take MRSA' structure functions |]12|; other modern parametrizations give very 
similar results. We use the leading order expression for a^, with Aqcd = 0.2 GeV, and take 
the (average) partonic px as factorization and renormalization scale. We use exact leading 
order matrix elements to compute the backgrounds from 2 — > 4 processes. These have been 
computed in ref.|T^ for the 4-jet final state, in ref. fl^ for the jjj'j final state, in ref.|T^ for 
JJ77 production, and in ref. [0 for jjee production. 
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In order to approximately mimic the acceptance of the CDF and DO detectors, we require 
all jets to have rapidity |?/jet| < 3.5, while we require Ive^l ^ 2.5 for electrons and photons. We 

also require the isolation cut ^Rij = \J {yi — Ujf' + (0i — > 0.7 for all combinations ij of 
final state particles. We generally find that the 4 — > 4 signal decreases more quickly than the 
2^4 background when the (transverse) momentum of the outgoing particles is increased. 
The reason is that the signal cross section contains four factors of parton densities, while the 
background only has two. We therefore try to keep the minimal acceptable px as small as 
possible, subject to the constraint that the event can still be triggered on. Specifically, we 
chose 

t) for 4-jet: Pt(Ji, J2) > 20 GeV, Pt(J3, J4) > 10 GeV 
ii) for iii-i: prilJi) > 15 GeV, ]5T(j2, ja) > 10 GeV; 
in) for jj77: Pt(7i, 72, ji, J2) > 10 GeV; 
iv) for jjee: pT(ei,e2) > 15 GeV, PT(ji,j2) > 10 GeV. 

The signal and background cross sections with only these basic acceptance cuts included 
are listed in column 2 of Table 1 for the 4-jet and jjj'j final states, and Table 2 for the j 777 
and jjee final states. We see that without further cuts, 4^4 processes only contribute 
between 7% (4-jet) and 18% (jjee), so additional cuts are clearly needed to extract the 
signal. As expected from our previous discussion, the signal to background ratio is worst for 
the 4-jet final state. 

As mentioned earlier, in 4 4 processes two pairs of particles are produced with equal 
and opposite transverse momenta, pr(l) = — Pt(2) andpT(3) = — pr(4). However, additional 
radiation can change the kinematics significantly, and the finite resolution of real detectors 
means that we can require momenta to be equal only within the experimental uncertainty. 

In the presence of initial or final state radiation the transverse momenta within a pair 
no longer balance exactly even if the resolution was perfect. We include this effect only 
for the signal, since in the background the final state particles in any case only pair up 
"accidentally"; we do therefore not expect large effects on the backgrounds. We randomly 
generate transverse "kicks" for each of the 2 — 2 processes in the signal. We assume that 
the direction of the kick is not correlated with the plane of the hard scattering. The absolute 
values qt of these additional transverse momenta are generated according to the distribution 

fiqr) oc exp [- {qo/qrf''^] /ql, (3) 



with < < qT.ma.x- This function describes the transverse momentum distribution 



of W bosons produced at ^/s = 1.8 TeV quite well, with go = 9 GeV. We adopt this choice 
of go for the jjee final state, which is dominated by the production of real Z bosons, but 
use the smaller value go = 4.5 GeV for the other final states, which are characterized by a 
smaller momentum scale. Finally, we take gT.max = 8 GeV as our default value; this assumes 
that one can reliably veto against jets with transverse momentum exceeding this value. 

We simulate finite energy resolutions by fluctuating the energies of all outgoing particles 
(keeping the 4-vectors light-like), using Gaussian smearing functions. The width of the 
Gaussian is given by 

5{E) =a-^/E®b-E, (4) 
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where © stands for addition in quadrature and E is in GeV. We take 

ajet = 0.80, 6jct = 0.05, a^,^ = 0.20, 6e,7 = 0.01, (5) 

which roughly corresponds to the performance of the CDF detector. We do not fluctuate 
the directions of the outgoing particles in this step. These are, however, affected by the 
transverse "kicks" mentioned earlier. For this reason, and in order to allow for an error in 
the determination of jet axes, we apply a relatively mild cut on the azimuthal opening angle 
of each pair: 

cos(0i - 0j) < -0.9. (6) 

This allows an opening angle as small as 154°. As emphasized earlier, in 4 — > 4 processes, 
the members of a pair should also have equal absolute values of px- As our final cut, we 
therefore require 

I \PTi^ I - IptU) 1 1 < c,,^J5WTi^)]+5VTiM, (7) 

with 5(|pr|) = ci ■ \J\pt\ © b ■ \pt\ as in eqs.@ and (|^). 

Our results for signal and background with these additional cuts included are summarized 
in the Tables. For the 4-jet and jjj'y final states (Table 1) we always take cu = cu = c, 
but we occasionally allow Cee^-yy > Cjj in the jjee and jj'j'j final states. The reason is that 
the cut is much more severe for e~^e~ and 77 pairs than for jet pairs, due to the better 
resolution of electromagnetic calorimeters, see eq.(^. Inclusion of the transverse "kick" 
therefore leads to a significant loss of signal if we take Cee,77 = 1. Although the stronger 
cut still gives a slightly better signal to noise ratio, given the limited available event sample 
employing a looser cut might give a statistically more significant signal. We do not attempt 
to quantify this statement here, since we have not included any reconstruction efficiencies in 
our calculation. Finally, in the last three columns of Table 1 we increase the cut on ARij 
from 0.7 to 1.2. This enhances the signal to background ratio by about 20 to 25%. 

Switching on energy smearing and transverse momentum kicks, and imposing the cuts 
and (0) with c = 5, reduces the signal by typically a factor of 2. This reduction is 
almost entirely due to the energy smearing. Ignoring the transverse kicks for the moment, 
in the signal both members of a pair have equal \pt\- If it falls below the cut-off value, both 
energies have to fluctuate upwards for the event to be accepted. In contrast, the downwards 
fluctuation of one energy can be sufficient to remove an event from the sample. The reduction 
is smaller for jjee production since most electrons have typically px — Mz/2, well above the 
lower limit. Fortunately the background is reduced even more in this step, by a factor of 4 
for 4-jet and jjj'j and 9 for jjee and jj77 final states, mainly due to the cut (H). Making 
the cut (|^ stricter, i.e. decreasing c, only slightly enhances the signal to background ratio 
in Table 1. This is partly due to the transverse kicks. Without them, the jjjj signal for 
c = 1.0 would be about 50% larger. This indicates that restricting additional jet activity as 
much as possible is quite important. 

Although in Table 1 the optimized S/B ratios are only about 0.23 for 4-jet and 0.31 for 
jjj'j, the signals are statistically quite significant; recall that the CDF and DO experiments 
together have accumulated about 200 pb~^ of data. We do not attempt to further optimize 
the S/B ratio for these two processes because we do not trust our parton level analysis, with 
a simplified treatment of finite detector resolutions and the effect of parton showering, suf- 
ficiently to extrapolate into the tails of distributions. Nevertheless, given the normalization 
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uncertainties of leading order QCD predictions, one will have to study the shapes of various 
distributions, such as the opening angle cos0ij, Ai?jj and pr balancing etc. in order to con- 
vince oneself that a signal is indeed present. Clearly the S/B ratio is much more favorable 
for the jjee and jj77 final states (Table 2). For these final states reducing c from its starting 
point c = 5 does increase this ratio significantly. Recall that for a fixed value of c the cut 
is much more restrictive for e~^e~ and 77 pairs than for jj pairs; this reduces the background 
more than the signal. On the other hand, this also has the effect that after imposing the cut 
(|^ with Cjj = c-y^ = 1, the size of the jj77 signal depends quite sensitively on the treatment 
of the transverse kick. Had we used go = 9 GeV in eq.(|), as appropriate for W production, 
the signal would have been reduced by a factor of about 0.7, while without any transverse 
kick it would have been larger by a factor 1.6. Clearly this uncertainty can be reduced by 
using the actual measured pt distribution of 77 pairs produced at the Tevatron. Fortunately 
the jjee signal is less sensitive to the "kick", since the electrons are usually so hard that 
adding or subtracting a few GeV does not matter very much. This final state therefore offers 
our most promising and robust signal. 

In summary, we have studied four different final states with a view of establishing an 
unambiguous signal for multiple partonic interactions in pp collisions at the Tevatron. The 
4-jet and jjj'~f final states offer very large event samples, but with a S/B ratio about 0.2 - 
0.3. One must study the shapes of various kinematical distributions for confirmation of the 
existence of the signal, as was indeed done by the CDF collaboration in their study of the 
4-jet final state 0. The situation is much more favorable for the jj'j'j and, especially, jjee 
final states; in the latter case one can increase the event sample by including muon pairs as 
well. Although even in these channels the signal to noise ratio is less favorable than what 



we found for four-jet production in 77 collisions a clear signal should be visible already 
in the present data sample. 

Once a signal is found, it would be important to establish if the normalization (Tq in 
eq.([^) is indeed the same for different processes, and independent of the Bjorken— a; range 
probed, as assumed in minijet models. Further, it would be very interesting to reduce the pt 
cut for at least some of the jets as much as possible, so that one can get closer to the actual 
minijet region. This could greatly enhance our understanding of "minimum bias" physics, 
and give us some confidence that we can trust extrapolations to LHC energies, where the 
understanding of overlapping minimum bias events becomes a crucial issue in the assessment 
of the viability of various "new physics" signals. Finally, such studies might shed new light 
on the thirty-year old problem of the rising total hadronic cross sections. 
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Table 1: Signal and background cross sections, as well as their ratios (S/B), for 4-jet 
production (in nb) and jjj'y production (in pb) at the Tevatron. In the first column only 
the basic acceptance cuts on the transverse momenta, rapidities and on ARij have been 
applied. In the second column we in addition apply the cuts (|D and (0), with c = 5. In 
the last three columns we sharpen the AR cut to ARij > 1.2, and gradually reduce c as 
indicated. Note that the "basic" cross sections have been computed ignoring finite energy 
resolution and transverse "kicks"; these effects have been included in the other columns, as 
described in the text. 







ARy > 0.7, 


ARij > 1.2 




basic 


c = 5 


c = 5 


c = 2 


c = 1 


^(4j)(S) 


266 


131 


91 


87 


57 


a(4j)(B) 


3,990 


878 


485 


442 


246 


S/B 


0.067 


0.15 


0.19 


0.20 


0.23 


f^(jJJ7)(S) 


515 


265 


169 


158 


97 


(^(jJJ7)(B) 


5,370 


1,310 


611 


571 


311 


S/B 


0.096 


0.20 


0.28 


0.28 


0.31 



Table 2: Signal and background cross sections in pb, as well as their ratios, for jj'yy 
production and jje^e^ production at the Tevatron. The notation is as in Table 1, except 
that we use the basic isolation cut ARij > 0.7 everywhere, and allow different values for 
Cl = Cjj and C2 = Cge or c^«,. 





basic 


Cl = C2 = 5 


Cl = C2 = 2 


Cl = 1, C2 = 2 


Cl = C2 = 1 


o-(jJ77)(S) 


1.86 


0.96 


0.71 


0.59 


0.37 


(^(jJ77)(B) 


20.8 


2.34 


1.16 


0.94 


0.52 


S/B 


0.089 


0.41 


0.61 


0.63 


0.71 


o"(jjee)(S) 


3.45 


2.01 


1.42 


1.07 


0.62 


t^(jjee)(B) 


19.0 


1.94 


1.00 


0.70 


0.37 


S/B 


0.18 


1.04 


1.42 


1.53 


1.68 
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